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turbulence  modeling  techniques  should  be  re-examined. 
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I.  INTRODUCTION 


The  turbulent  wake  produced  by  a body  in  a fluid  stream  has  Iona 
been  a subject  of  interest  in  fluid  mechanics.  In  more  recent  years, 
with  the  increased  concern  for  maintaining  the  natural  ecological 
system  of  checks  and  balances,  this  topic  has  acquired  added  impact. 
With  a thorough  understanding  of  wake  development  and  subsequent  dissi- 
pation, the  problem  of  dispersion  and  control  of  dispersion  of  pollu- 
tants from  aircraft  and  water  vehicle  propulsion  units  can  better  be 
approached.  Experimental  studies  of  the  turbulent  wake  are  also  of 
interest  for  testing  theories  on  turbulent  modelling. 

Unfortunately,  as  is  evident  from  Table  I,  most  of  the  previous 
experimental  studies  have  dealt  with  relatively  simple  bodies  such 
as  circular  disks,  ellipsoids,  spheroids,  and  some  slender  bodies. 
Generally  these  studies  have  been  restricted  to  the  unpropelled  or 
"drag-body"  situation.  In  fact,  most  estimates  of  the  wake  develop- 
ment behind  self-propelled  bodies  have  been  based  on  the  experiments 
of  Naudascher  [1]  in  which  the  momentumless  condition  was  established 
with  a circular  disk  and  a central  jet  providing  the  necessary  thrust 
to  balance  the  drag  of  the  disk.  Later,  Gran  [2]  investigated  the 

momentumless  wake  behind  a self-propelled  Rankine  ovoid  at  Re^  = 

4 

6 x 10  . Until  recently  the  above  were  the  only  two  self-propelled 
configurations  which  had  been  examined  and  the  utility  of  these  re- 
sults are  limited  because  of  impractical  body  shapes  or  low  Reynolds 
numbers,  or  both.  More  extensive  work  was  published  in  1974  by 
Swanson  [3],  et.al.,  and  Chieng  [4 | , et.  al.,  at  Re0  * 6.18  x 105. 
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TABLE 


Summary  of  Subsonic,  Axisymmetric  Turbulent  Wake  Experiments 
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These  studies  provide  a systematic  comparison  of  the  turbulent  wakes 
behind  slender  bodies  with  identical  forebody  shape  by  varying  through 
the  following  cases:  1)  pure  drag  body,  2)  self-propelled  by  axial 

fluid  injection,  3)  self-propelled  with  a well  designed  propellor. 

Daffan  [5]  extended  the  work  of  Swanson  and  Chieng  on  the  propellor- 
driven  configuration  by  studying  the  effect  on  the  wake  of  a sail -type 
appendage  to  the  formerly  axisymmetric  body  and  also  the  effect  of 
nonzero  pitch  angles  (0°  to  -2°).  Schetz  [6],  et.al.,  studied  the 
effect  of  replacing  the  single  propellor  by  an  equivalent  set  of  side- 
by-side  counter  rotating  propellors.  Additionally,  Schetz  conducted 
measurements  of  the  aero-hydrodynamic  forces  and  moments  on  propel lor- 
driven  slender  bodies,  and  boundary  layer  measurements  on  the  body 
itself. 

It  is  evident  that  the  work  reported  in  Refs.  No.  [3],  [4],  [5], 
and  [6]  has  contributed  greatly  to  the  better  understanding  of  the 
phenomena  associated  with  the  turbulent  wakes  produced  by  self-propelled 
bodies  of  practical  configurations.  The  present  investigation  was  de- 
signed to  study  the  development  of  the  mean  wake  flow  using  the  fluid 
temperature  as  a tracer  and  to  document  the  turbulent  fluctuations 
including  those  of  temperature.  The  stratification  of  the  temperature 
field  remains  sufficiently  small  that  buoyancy  effects  are  negligible. 

The  testing  was  conducted  in  the  VPI  6'  x 6‘,  subsonic  wind  tunnel 
5 

at  Rey  ' 2.04  x 10  . The  approaching  mean  velocity  and  temperature 
fields  were  measured  with  a pitot-static  tube  and  a thermocouple. 

Mean  boundary  layer  measurements  were  made  at  several  stations  along 
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the  body  with  a "rake"  of  small  pitot  tubes  for  the  velocity  field 
and  a small  thermocouple  for  the  temperature  field.  The  mean  velocity 
and  temperature  fields  in  the  near  wake  were  measured  with  a yawhead 
probe  and  a thermocouple  respectively.  From  the  yawhead  probe  the 
flow  angularity  could  also  be  determined.  Measurements  of  axial  and 
radial  turbulence  intensities,  radial  shear  stress,  and  temperature 
fluctuations  were  taken  using  cross-wire  and  straight  hot-wire  probes. 

In  order  to  achieve  a temperature  stratified  flow,  an  air  injection 
and  heating  system  was  designed  and  built.  The  primary  design  cri- 
terion of  the  heating/ injection  arrangement  was  to  produce  in  the  wind 
tunnel  main  stream  a slug  of  relatively  high  temperature  air  while  at 
the  same  time  disturbing  the  normal  dynamic  pressure  field  as  little 
as  possible. 

This  report  is  arranged  into  several  major  sections.  First,  the 
apparatus,  models,  and  experimental  methods  are  described  in  detail. 

The  results  are  presented  in  separate  sections,  starting  with  the 
boundary  layer  data.  The  next  subsections  in  the  RESULTS  section  deal 
with  the  mean  wake  flow  and  finally  the  turbulence  measurements  in  the 
near  wake  are  reported.  The  results  are  presented  in  graphical  and  then 
tabulated  form  in  Appendices  at  the  rear  of  this  report. 


II.  EXPERIMENTAL  APPARATUS  AND  INSTRUMENTATION 


Wind  Tunnel  Facility 

All  tests  were  conducted  in  the  VPI  & SU  6'  x 6'  subsonic  sta- 
bility tunnel  (Fig.  No. 1)  at  a dynamic  pressure  of  1.0  in.  of  water 

(approximately  68  ft/sec)  corresponding  to  a Re^  based  on  diameter  of 
5 

2.04  x 10  . This  facility  is  a continuous, closed  jet, single  return, 
wind  tunnel.  The  test  section  is  28  ft.  long  which  is  sufficient, 
with  the  model  located  near  the  center  of  the  section,  to  allow 
measurements  in  the  near  wake.  The  airstream  has  a low  turbulence 
intensity  factor  of  1.08.  To  control  the  tunnel  velocity  a Barocel 
Electronic  Manometer  is  used  to  read  dynamic  pressure  from  a pitot- 
static  tube  mounted  on  the  test  section  wall  out  of  the  wake  of  the 
model.  Static  free-stream  pressure  is  measured  by  a Validyne  digital 
barometer  Model  DB99. 

Sel f-Propel led  Model 

The  model  considered  in  this  investigation  has  an  overall  length 
of  72  in.  and  a maximum  diameter  of  6 in.,  giving  a fineness  ratio  of 
12:1.  The  model  was  strut  mounted  from  the  floor  of  the  wind  tunnel 
in  the  center  of  the  test  section  (Fig. No.  2).  The  forebody  is  para- 
bolic and  machined  from  laminated  layers  of  plexiglass.  The  center- 
body  is  an  aluminum  cylindrical  tube.  The  tail  body  is  also  plexiglass 
machined  into  an  ogival  shape.  The  overall  dimensions  of  the  body 
are  shown  in  Fig.  No.  (3). 

The  model  used  a 2.75  hp.  DC  motor  to  direct  drive  the  single 
propel  lor  shaft  which  extended  thru  the  center  of  the  tail  section 
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( 

(Fig.  No.  4) . A 6 in.  diameter,  4 in.  pitch  three-bladed  model  air- 
plane propellor.  Tornado  Model,  manufactured  by  Griggs  Bros,  was 
used  to  provide  sufficient  thrust  to  obtain  the  self-propelled 
condition. 

To  avoid  overheating  in  long  duration  experiments,  the  motor 
was  enclosed  by  copper  tubing  to  allow  for  water  cooling.  The  water 
cooling  was  not  used  for  the  short  duration  test  necessary  to  de- 
termine the  required  propellor  speed  to  produce  the  momentumless 
condition  since  the  motion  of  the  cooling  water  and  the  connecting 
tubing  could  affect  this  determination.  The  "sail-drag"  is  included 
in  determining  this  condition.  In  order  to  measure  the  drag  forces 
on  the  body  and  subsequently  measure  the  thrust  necessary  to  exactly 
counteract  these  forces,  a strain  gage  balance  (Fig.  No.  5)  was 
internally  mounted  in  the  model.  To  power  the  DC  motor,  a Sorenson 
150v,  15a  power  supply  was  used. 

Ai r Heating  Ajjpar atus 

The  stratification  of  the  approaching  airstream  was  induced  by 
injecting  heated  air  at  the  freestream  velocity  upstream  of  the  self- 
propelled  model.  A schematic  of  the  test  section  arrangement  includ- 
ing the  injection  and  self-propelled  model  is  shown  in  Fig.  No.  (6). 
For  this  purpose  an  injector  was  constructed  with  the  overall  dimen- 
sions shown  in  Fig.  No.  (7).  The  interior  of  the  injector  is  shown 
in  Fig.  No.  (8).  Air  enters  the  injector  through  the  smaller  1 in. 
galvanized  pipe  and  exits  through  1/8  in.  holes  which  are  distributed 
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uniformly  on  the  centerline  along  the  span.  Surrounding  this  pipe  is 
a larger  2 in.  galvanized  pipe  into  which  two  rows  of  1/2  in.  holes 
are  evenly  distributed  along  the  span.  These  holes  are  offset  above 
and  below  the  centerline  by  45°.  Downstream  of  the  pipe  assembly, 
there  is  a .035  in.  stainless  steel  plate  to  help  guide  the  exiting  air 
and  begin  the  flow  straightening  process.  Immediately  downstream  of 
this  plate,  the  injector  is  filled  with  1/8  in.  aluminum  tubes  of  1.25 
in.  length.  These  are  held  in  place  by  a small  box  with  sides  and  top 
constructed  of  .035  in.  stainless  steel.  The  upstream  and  downstream 
ends  of  the  box  are  covered  with  .011  in.  copper  mesh.  The  mounting 
of  the  injector  may  be  seen  in  Fig.  No.  (2).  The  basic  ideas  for 
the  construction  of  the  two-dimensional  injector  were  obtained  from 
Hokenson  [9], 

The  injected  air  was  supplied  by  four  Ingersol 1 -Rand  Type  90 
compressors  ooerating  continuously.  The  compressed  air  fed  into  a 
Grove  dome  pressure  regulator  which  controlled  the  mass  flow  rate.  The 
heating  was  done  electrically  by  6 Plasmatron  PS-20  transformers.  The 
heating  mechanism  was  as  shown  in  Fig.  No.  (9).  After  leaving  the 
regulator  the  air  was  guided  into  two  lengths  of  Inconel  601  1-1/4 
in.  diam.,  1/4  in.  wall  pipe.  The  Plasmatrons  were  connected  to  the 
upstream  and  downstream  ends  of  the  pipes  thus  forming  a parallel 
electric  circuit.  To  minimize  the  heat  loss  to  the  atmosphere,  the 
entire  piping  system  from  the  upstream  end  of  the  heating  circuit  to 
the  injector  was  covered  by  1-1/4  in.  Kay-Lo  Owens-Corning  insulation. 


L . . . .. - A 


Wake  Probes 


The  mean  flow  pressure  measurements  in  the  wake  were  taken  using 
a three-dimensional  yawhead  probe  constructed  by  United  Sensor  (Fig.  No. 
10).  The  output  from  the  pressure  probe  was  measured  on  a Barocel 
electronic  manometer  and  displayed  on  a Doris  digital  voltmeter.  To 
read  all  five  pressure  points  for  the  yawhead  probe  on  the  manometer, 
a Scanivalve  Model  W0601/IP-RT  fluid  switch  wafer  with  solenoid  drive 
was  used.  The  mean  temperature  field  was  measured  by  means  of  an 
Omega  SCASS-062G-6  thermocouple.  The  output  was  fed  into  a Consoli- 
dated Controls  Type  K digital  indication  voltmeter  and  read  directly 
in  °F. 

The  axial  and  radial  trubulence  and  radial  shear  stress  data  were 
obtained  using  a cross-wire  probe  (TSI  Model  1241T1.5).  The  probe 
was  operated  at  an  overheat  ratio  of  1.8  using  a constant  temperature 
anemometer  module  (TSI  Model  1050)  and  a power  supply  and  monitor 
(TSI  Model  1051-6).  Readings  were  taken  from  a DISA  55D35  RMS  meter. 

A Model  1015C  correlator  and  TSI  Model  1076  voltmeter  were  also  used. 
Temperature  fluctuations  in  the  wake  were  measured  using  a straight  hot 
wire  (TSI  Model  1210).  This  probe  was  operated  at  a low  overheat  ratio 
of  1.05  in  order  to  more  accurately  measure  the  temperature  fluctuations. 

Traverse  Mount 

The  traverse  mount  had  a vertical  movement  of  4 ft.  and  a hori- 
zontal movement  of  11  in.  Both  vertical  and  horizontal  movement  was 
accomplished  with  variable  speed  motors  which  were  remote-controlled 
from  outside  the  test  section.  The  probe  position  was  monitored  using 
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a series  of  potentiometers. 


Boundary  Layer  Probes 

The  pitot  rake  (Fig. No.  11 ) was  developed  by  Nerney  [7]  and  was 
approximately  1.5  in.  high  with  one  static  port  and  24  total  pressure 
tubes  at  increasing  distances  up  the  rake.  This  allowed  for  more  data 
points  in  the  wall  region.  The  temperature  field  near  the  body  was 
measured  by  means  of  an  Omega  SCASS-062G-6  thermocouple  suspended 
vertically  from  the  traverse  mounts  . Ih^ thermocouple  could  then  be 
withdrawn  to  any  desired  distance  from  the  body. 


III.  EXPERIMENTAL  METHODS 


Producing  a Nearly  Uniform  Velocity,  Temperature-Stratified  Flow 
Initially,  the  primary  consideration  was  to  define  the  control 
settings  which  would  allow  a nearly  uniform  velocity  distribution 
downstream  of  the  injector  while  maintaining  (T^  - T)  as  large  as 
possible.  This  was  accomplished  in  the  following  sequence.  The 
Plasmatrons  were  activated  at  the  highest  setting  which  they  could 
safely  support.  The  compressors  were  then  activated  and  allowed  to 
run  continuously  forcing  the  air  thru  the  heating  system  and  thru  the 
injector.  The  tunnel  was  operated  at  the  test  condition  of  1.0  in. 

H^O  dynamic  pressure  (approximately  68  ft/sec).  The  injector  exit 
temperature  was  monitored  by  a thermocouple  located  .5  in.  from  the 
exit.  After  about  20  minutes,  this  temperature  would  stabilize,  in- 
dicating the  steady-state  condition  had  been  reached.  At  this  time 
temperature  and  pressure  scans  were  taken  at  a station  16  in.  upstream 
of  the  submarine  body.  The  velocity  defect  or  overshoot  was  observed 
and  appropriate  adjustments  were  then  made  to  the  pressure  regulator 
ahead  of  the  heating  system.  By  trial  and  error,  the  proper  control 
settings  were  determined  to  be  as  follows:  The  compressors  were  oper- 

ated continuously  between  120-130  psig.  The  regulator  was  set  to  main- 
tain a 60  psig.  exit  pressure.  The  large  pressure  difference  between 
the  incoming  flow  to  the  regulator  and  the  exiting  flow  was  necessary 
to  help  facilitate  a constant  flow  through  the  regulator.  The  plasma- 
trons were  operated  at  1740a  and  20v.  These  settings  resulted  in  nearly 
uniform  velocity  profiles  and  substantial  temperature  nonuniformity  at 
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the  test  station  ahead  of  the  submarine  body.  They  corresponded  to  a 
U.j/Ux  of  approximately  1.3  and  a of  approximately  150°F. 

Typical  velocity  and  temperature  profiles  at  the  station  16  in.  up- 
stream of  the  submarine  body  are  shown  in  Fig.  No.  (12).  Note  the 
greatly  exaggerated  scale  for  velocity  variations.  Precise  lateral 
uniformity  from  the  injector  proved  to  be  difficult  to  obtain.  The 
basic  magnitudes  of  the  velocity  and  temperature  profiles  were  later- 
ally fairly  uniform,  but  the  local  maximums  were  often  at  different 
locations  along  the  Y-axis.  This  "undesigned"  characteristic  of  the 
injector  proved  to  be  helpful,  however,  in  tracing  the  wake  flow. 

Producing  the  Self-Propelled  Condition 

To  determine  a self-propelled  configuration, a strain  gage  balance 

was  internally  mounted  in  the  model  (Fig. No. 5).  With  the  model  mounted 

in  the  tunnel,  potentiometers  were  used  to  zero  the  output  from  the 

balance  due  to  the  weight  of  the  model.  The  injector  was  operated,  and 

the  free-stream  velocity  was  increased  to  1.0  in.  H^O  while  the  axial 

force  was  monitored.  The  propel  lor  rpm  was  then  increased  until  the 

axial  force  reading  was  zero.  It  was  determined  that  a propel  lor  rpm 

of  12987  was  needed  to  achieve  the  self-propelled  status.  This  corre- 

/U  \ 

sponds  to  advance-diameter  ratio!  of  0.628. 

Boundary  Layer  Measurements 

Measurements  were  made  with  the  pitot  rake  and  the  thermocouple  at 
Z'/D  = 2.67,  4.67,  6.67,  8.67,  and  10.0.  The  data  were  reduced  using 
the  Bernoul  1 i equation,  the  local  pitot  pressure,  and  the  local  static 
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temperature  obtained  from  the  thermocouple.  The  local  static  pressure 
was  determined  by  lowering  a small  pitot-static  tube  to  a point  0.25 
in.  from  the  body  and  not  from  the  static  tube  on  the  rake.  The 
boundary  layer  measurements  were  made  on  the  side  away  from  the  sail. 

The  coordinate  system  used  for  these  measurements  is  illustrated  in 
Fig  . No.  (13). 

Wake  M easurements 

The  mean  flow  velocity  as  well  as  the  flow  angularity  were  determi- 
ned using  a five-port  yawhead  probe.  This  probe  was  calibrated  over  a 
range  of  flow  pitch  and  yaw  angles  between  -30°  and  +30°  by  Daffan  [5]. 

It  was  decided  that  an  adeguate  description  of  the  wake  could  be 
achieved  with  three  traverses  which  are  shown  in  Fig.  No.  (14).  The 
three  traverses  were  made  at  downstream  stations  of  Z/D  = 0.33,1.0,  2.0, 
3.0  and  4.0.  The  coordinate  system  for  the  wake  study  is  illustrated 
in  Figs . No.  (13)  and  (15) . 

To  obtain  the  mean  velocity  components  from  the  pressure  measure- 
ments the  following  equation  for  total  pressure  was  used: 

PQ  = Ps  + 1/2PV2  + 1 /2pv 1 2 

where  V is  the  mean  total  velocity  and  vT  is  the  sum  of  the  mean 

2 2 2 

fluctuations.  Since  V > > v‘  the  1 / 2pv ' term  was  neglected.  The 
yawhead  probe  gave  direct  readings  of  P^  - P^,  where  P^  is  the  mean 
total  pressure  and  Pc  is  the  free-stream  static  pressure.  Therefore 


= Ps  + l/2pv2 


since 
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the  mean  total  velocity  can  be  calculated  from 

V, 


'TOTAL 


J 2 HP,  - ps,)~ 


The  mean  axial  velocity  can  be  obtained  by 


U, 


Vcos(FP)cos(FY) 


where  FP  is  the  flow  pitch  angle  and  FY  is  the  flow  yaw  angle.  Similar- 
ly the  other  velocity  components  are  given  by, 

Ux  = Vcos(FP)sin(FY) 

Uy  = Vsin(FP) 

To  obtain  the  local  static  pressure  from  the  yawhead  probe  a 
method  described  by  Winternitz  [8]  was  used.  Since  this  method  is  reliable 
only  for  substantial  flow  angularity  in  one  direction,  calibration  curves 
were  used  for  both  pitch  and  yaw.  The  flow  yaw  angle  FY  is  found  from 
the  calibration  probe  constant 

K - 

FY  ' Pt  - P4 

for  each  station.  Defining  the  hole  constant  as 

P,  - P, 


Ki (FY) 


J S 

q 


where  P.  is  the  hole  pressure  and  q is  the  dynamic  pressure,  Winternitz 


showed  that 


q = 


<P,  - PS->  - <P2  - PS-> 
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or 


(RS  ' P$J  (P2  ‘ PsJ  ‘ qK2 


Therefore 


<ps  - ps  > 


K1  (P2  - PsJ  - K2  (P1  - PS-> 
K,  - K2 


Similarly  for  the  pitched  case  where  the  yaw  angles  are  small 


TP 


P4  - P5 
P1  - P2 


and. 


„ „ K1  <PS  - PS«>  - K5  (p,  - PS.) 

S " s“  K1  ' K5 

From  these  relations,  the  local  static  pressure  could  be  determined 
when  there  were  substantial  pitch  angles  with  small  yaw  angles  or  for 
large  yaw  angles  with  small  pitch  angles.  Figure  (16)  is  an  illustra- 
tion of  the  sign  conventions  used  for  measuring  pitch  and  yaw. 

Mean  temperature  rreasurements  were  taken  in  the  wake  by  using  a 
thermocouple  and  making  the  identical  traverses  at  each  downstream 
location  behind  the  propel  lor.  Temperature  was  read  directly  in  °F  on 
a digital  vol tmeter . 

The  axial  and  radial  turbulence  intensities  and  radial  shear  stress 
were  measured  at  Z/D  =0.33  and  Z/D  = 4.0.  These  traverses  were  made  at 
each  station  as  was  done  in  the  mean  flow  study.  The  cross-wire  probe 
was  operated  at  a high  overheat  ratio  of  1 .8  in  order  to  decrease  the 
temperature  sensitivity.  Initially  the  probe  was  calibrated  by 
orienting  one  of  the  sensors  parallel  to  the  flow  direction  and 
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subjecting  the  probe  to  two  different  flow  conditions  (motionless  air 
and  air  moving  at  98.4  ft/sec).  The  probe  was  then  rotated  90°  to 
align  the  other  sensor  and  the  process  was  repeated.  Temperature  was 
recorded  for  each  situation. 

Figure  (17)  shows  the  two-sensor  orientation.  With  this  arrange- 
ment the  turbulence  in  two  perpendicular  directions  and  the  directional 
correlation  (Reynolds  stress)  can  be  measured.  Assuming  that  the 
effective  cooling  velocity  follows  a sine  law  and  that  the  two  sensors 
are  45°  to  the  mean  flow,  their  output  voltages  are  related  by: 

A = SA  (U  + V) 

B = SB  (U  - V) 

where  S = proportionality  constant,  A and  B are  output  voltages.  If 
the  two  signals  are  adjusted  so  that  = Sg  = S (during  calibration) 
then: 


Vu'2  = 

X 1 

2S 

1 1 (A  + B)L 

. 

h i 

J (A  - B)2 

uV  = 

f2 

(A2  - B2) 

where  f appears  as  a temperature  compensating  factor  which  corrects  for 
temperature  variations.  For  the  measurements  of  this  experiment  the 
calibration  constants  were  determined  as  follows: 

S = .0718 

f = J 396.8  - T 
J 347.8 
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where  T = °F. 

Temperature  fluctuation  measurements  were  made  using  a straight- 
wire  probe  at  a low  overheat  ratio  of  1.05  to  minimize  velocity  sensi- 
tivity and  maximize  temperature  sensitivity.  The  probe  was  calibrated 
by  immersing  the  hot  wire  into  first  an  ice  bath  and  second  into  boiling 
water.  The  output  was  adjusted  so  that  the  measured  voltage  was  O.Ov 
at  32°F  and  lOv  at  212°F.  Output  voltage  varies  linearly  with  temper- 
ature between  these  extremes.  The  probe  was  then  mounted  on  the 
traverse  and  the  wake  traverses  were  repeated  at  Z/D  =0.33  and  Z/D  = 4.0. 


IV.  RESULTS 


The  experimental  results  obtained  during  this  research  are  presented 
in  Figs.  No.  (18)  - (49)  in  the  form  of  velocity  and  temperature  pro- 
files for  each  station  examined  along  the  body.  Figures  (20)  - (34)  are 
graphs  of  the  distribution  of  flow  angle,  mean  static  pressure,  and 
axial  velocity  for  each  station  examined  in  the  wake.  Figures  (35)  - 
(39)  are  graphs  of  the  mean  temperature  distribution  at  each  of  these 
stations.  Figures  (40)  - (49)  are  graphs  of  the  axial  and  radial  tur- 
bulence intensity,  radial  shear  stress,  and  temperature  fluctuation  dis- 
tributions at  Z/D  = 0.33  and  Z/D  = 4.0.  The  results  are  also  presented 
in  tabular  form  in  Tables  A-l  thru  A-14. 

The  turbulence  data  were  not  corrected  for  the  effect  of  flow  angu- 
larity and  the  mean  flow  data  were  not  corrected  for  the  effect  of  tur- 
bulent fluctuations.  The  latter  approximation  is  valid  to  a high  degree 
of  accuracy  due  to  the  relatively  small  turbulence  levels  in  most  of  the 
wake  region.  Jakubowski  (23)  has  found  in  recent  tests  that  the  error 
induced  by  not  correcting  the  data  for  flow  angle  variations  is  less  than 
7%  in  general  for  the  flow  angles  encountered  in  this  research.  Follow- 
ing is  a discussion  of  the  results  of  the  current  research. 

Body  Boundary  Layer  Observations 

The  boundary  layer  development  along  the  body  was  examined  in  order 
to  trace  the  temperature  distribution  in  the  longitudinal  direction  and 
subsequently  more  accurately  determine  the  effects  of  the  propellor  on 
the  near  wake  flow  field.  In  Fiq.  No.  (18)  are  plotted  the  boundary 
layer  velocity  data,  non-dimensional ized  by  U , for  Z 1 /D  = 2.67,  4.67, 

oo 
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6.67,  8.67,  and  10.0.  These  data  are  very  similar  to  those  first  re- 
ported in  Ref.  No.  (6).  They  show  a normal  boundary  layer  grow 
with  increasing  downstream  distance.  The  slope  of  the  profiles  near  the 
wall  tends  to  indicate  a developing  turbulent  boundary  layer.  An 
examination  of  the  boundary  layer  thickness  at  stations  Z‘/D  = 6.67, 

g 

8.67,  and  10.0  shows  that  s(x)  ~ x'  as  would  be  the  case  for  a tur- 
bulent boundary  layer.  Turbulence  measurements  were  not  made  in  the 
boundary  layer  due  to  lack  of  instrumentation  for  that  task.  These 
data  were  collected  with  the.  submarine  in  the  self-propelled  mode  and 
the  effects  of  the  propel  lor  on  the  development  of  the  boundary  layer 
were  not  studied  since  Schetz  (6)  has  determined  that  the  prop  has 
little  or  no  effect  until  Z'/D  > 11. 

Figure  No.  (19)  contains  the  temperature  distributions  near  the 
body  at  streanwise  stations  corresponding  to  those  of  Fig.  No.  (18). 
First,  it  should  be  noted  that  these  distributions  are  not  graphs  of  the 
thermal  boundary  layer  but  rather  that  the  thermal  boundary  layer  is 
contained  within  the  variable  temperature  region.  The  Prandtl  number, 
which  is  the  connecting  link  between  the  velocity  and  temperature  fields, 
is  near  unity  for  air  which  indicates  that  the  velocity  and  temperature 
boundary  layers  should  have  nearly  equivalent  dimensions.  The  tempera- 
ture varies  over  a normal  distance  of  nearly  double  that  of  the  velocity 
boundary  layer.  This  is  simply  a result  of  the  initial  injection 
geometry  and  subsequent  spreading  of  the  jet.  The  development  of  the 
thermal  boundary  layer  can  be  qualitatively  examined  by  focusing  atten- 
tion on  the  region  nearest  the  wall.  Note  that  these  profiles  resemble 
classical  insulated-wall  thermal  boundary  layers.  The  regions  wherein 


20 


dT/dy  ~ 0 can  be  seen  to  grow  very  distinctly  with  increasing  Z'/D  to 
the  point  that  the  size  of  the  region  at  Z'/D  = 10.0  is  nearly  twice 
that  at  station  Z'/D  = 4.67.  This  is  essentially  identical  to  the  be- 
havior of  the  velocity  profiles  of  Fig.  No.  (18).  The  thickening  of 
the  overall  layer  is  also  evident  in  Fig.  No.  (19),  but  the  growth 
due  to  the  proximity  of  the  body  and  the  normal  growth  rate  of  the  jet 
cannot  be  quantitatively  separated. 

Figure  No.  (19)  also  shows  the  decay  of  the  maximum  of  the  high 
temperature  layer  with  increasing  streamwise  distance.  This  energy  is 
lost  to  the  free  stream  and  manifests  itself  as  a general  overall  heat- 
ing of  the  entire  wind  tunnel  fluid.  Typically  during  a test  program 
lasting  30  minutes,  the  free  stream  temperature  would  increase  by  1°  to 

2°F.  The  fact  that  the  mass  flow  of  the  jet  was  very  small  compared  to 

m . 

the  total  mass  flow  ( — = .01)  and  the  presence  of  air  coolers  in  the 

tot 

wind  tunnel  circuit  kept  the  overall  heating  rate  very  low. 

Mean  Wake  Flow  Observations 

Figures  No.  (20)  - (34)  are  graphical  presentations  of  the  mean 
static  pressure,  mean  axial  velocity,  and  flow  angularity  distributions 
in  the  near  wake  at  stations  Z/D  = 0.33,  1.0,  2.0,  3.0,  and  4.0.  At 
each  station  results,  are  presented  along  three  distinct  traverses 
vertical,  horizontal,  and  e = 45°.  The  coordinate  systems  and  sign 
conventions  used  are  shown  in  Figs.  No.  (13)  - (16). 

In  general  from  an  overall  examination  of  these  data,  one  sees  the 
initial  momentum  excess  and  defect  regions  spread  in  extent  but  decrease 
in  magnitude  with  increasing  Z/D.  Similar  behavior  is  exhibited  by  the 
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propellor  swirl  and  static  pressure  variations.  Along  a vertical 
traverse  the  flow  is  pitched  very  slightly  and  the  yaw  angles  are 
larger.  The  converse  is  true  along  the  horizontal  traverse.  Maximum 
flow  angles  vary  from  approximately  10°  at  Z/D  = 0.33  to  about  5°  at 
Z/D  = 4.0.  Results  along  the  e = 45°  traverses  show  that  in  this  plane 
flow  pitch  and  flow  yaw  are  about  equally  present.  Some  small  asymmetry 
of  the  results  is  evident.  This  is  due  to  a combination  of  circumstances 
which  cause  the  flow  to  deviate  slightly  from  axisymetric  - the  presence 
of  the  sail,  the  vertical  and  lateral  gradients  in  the  injected  jet,  and 
the  small  errors  in  alignment  which  occur  in  positioning  the  model. 
Basically  the  behavior  of  the  flow  is  similar  to  that  reported  by  Daffan 
[5].  The  magnitude  of  the  momentum  excess  and  defect  and  the  other 
associated  properties  of  the  wake  structure  are  smaller  due  to  the  lower 
Reynolds  number  and  corresponding  lower  propellor  rpm  of  this  study.  One 
noticeable  departure  from  the  results  of  Ref.  No.  (5)  is  the  ap- 
pearance of  another  momentum  defect  region  in  the  near  wake.  At  every 
station  there  is  the  maximum  defect  in  the  wake  center.  Passing  the 
center  there  is  a region  of  momentum  excess,  characterized,  for  instance, 
by  the  axial  velocity  overshoot.  Continuing  into  the  outboard  portion 
of  the  wake  there  is  invariably  another  region  of  velocity  defect  which 
is  smaller  in  both  magnitude  and  extent.  Actually  this  behavior  is 
predicted  by  simple  propellor  theory,  and  is  discussed  in  Ref.  No. 

(24). 

Of  particular  interest  are  the  results  in  Figs.  No.  (35)  - (39). 
These  are  graphs  of  the  temperature  distribution  In  the  wake  at  the 
identical  stations  used  in  obtaining  the  mean  pressure  data.  At  each 
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station  the  three  traverses  were  made,  first  with  the  prop  not  rotating 
and  second  at  the  rpm  necessary  to  achieve  the  self-propelled  condition. 

In  the  former  circumstance,  the  prop  was  rotated  by  hand  before  the 
traverse  was  made  and  positioned  where  it  would  not  interfere  with  the 
flow. 

The  results  at  station  Z/D  = 0.33  ( Fig.  No.  (35))  are  typical  of 
each  station  if  viewed  singly.  The  initial  temperature  distribution  with- 
out the  influence  of  the  propeller  can  be  easily  seen.  The  vertical 
variation  is  quite  pronounced  with  peak  temperatures  of  approximately 
16°  - 20°  greater  than  T^.  Ideally  the  temperature  variation  along  the 
horizontal  traverse  should  be  zero,  however  variations  up  to  5°F  are 
observed.  This  is  largely  due  to  the  lateral  gradients  which  were  pre- 
sent initially  in  the  injected  jet  and, to  a lesser  degree, due  to  distor- 
tion produced  by  the  body.  The  local  temperature  maximums  at  each  lateral 
location  along  the  jet  injector  were  very  nearly  the  same,  however,  the 
vertical  location  of  these  maximums  (i.e.,  the  vertical  temperature  dis- 
tribution) would  vary  somewhat  laterally.  Rotating  the  propellor  has  the 
effect  of  mixing  the  flow  to  some  extent.  Regions  where  the  temperature 
was  higher  tend  to  decrease  in  temperature  and  virp-versa.  As  Z/D  in- 
creases, these  data  indicate  an  overall  rotation  of  the  flow.  This  is 
most  evident  from  the  vertical  and  horizontal  traverses  of  Figs.  No. 

(35)  - (39).  Before  activation  of  the  propel  lor,  the  highest  temperature 
levels  in  the  wake  of  the  body  are  located  along  the  X-axis  near  X = 1 .0 
in.  Examination  of  the  temperature  field  with  the  propellor  running 
shows  that  these  higher  levels  of  temperature  are  beginning  to  appear 
along  the  Y-axis  at  Y = -1.5  in.  further  downstream  at  Z/D  = 4.0. 
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Similarly  the  temperature  levels  along  the  x-axis  drop  considerably  to 
levels  that  are  consistent  with  those  that  existed  near  = 1.0  in.  before 
prop  activation. 

Consider  for  a moment  the  fluid  particles  that  existed  at  X =1.0 

o 

in.  and  Y = 0.00  in.  at  time  t = t . The  fluid  state  of  those  particles 
oo  r 

is  determined  by  the  data  illustrated  in  Figs.  No.  (21)  and  (35).  As 
has  been  pointed  out,  the  temperature  distributions  of  Figs.  No.  (35)- 
(39)  tend  to  indicate  that  they  have  moved  to  X 1 = 0.0  in  and  Y 1 ' -1.5 
in.  during  the  time  t-j  - tQ  necessary  for  them  to  negotiate  22.0  in. 
along  the  longitudinal  axis.  If  this  is  in  fact  the  case,  calculations 
based  on  the  mean  pressure  data  of  Figs.  (21)  - (35)  should  predict 

this  rotation  and  translation.  Unfortunately,  the  coarse  "grid  size" 

employed  for  the  data  collection  prevents  an  exact  analysis  of  the  path 
of  a particular  streamline,  but  by  taking  the  fluid  state  as  it  existed 
at  station  Z/D  = 0.33  and  assuming  that  all  gradients  of  the  fluid  pro- 
perties are  zero,  one  can  approximately  predict  the  location  of  these 
particles  at  time  t^  - tQ.  Carrying  out  the  indicated  calculations  the 
results  are, 

Uy  = -7.25  ft/sec 

' lfx  = =1 ‘.20  ft/sec  - - - - - 

t,  - t = 0.026  sec 
I o 

Y1  - Yo  = -1.5  in 

X.  - X = -0.382  in 
1 o 

The  results  of  the  approximate  calculations  indicate  that  these  particles 
would  move  in  pitch  more  than  is  observed  and  move  in  yaw  less  than  is 
observed.  The  direction  of  rotation  is,  of  course,  the  same  as  that 
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observed  (counter-clockwise  looking  upstream).  The  discrepancies  be- 
tween the  calculations  and  the  observed  results  are  due  to  the  assump- 
tion of  zero  gradients  in  the  flow  properties.  In  fact,  as  the  fluid 
particles  move  from  the  +X  axis  to  the  -Y  axis,  the  pitch  angles  become 
smaller  and  the  yaw  angles  become  larger.  This  would  result  in  a fluid 

element  moving  more  in  yaw  and  less  in  pitch  than  is  crudely  predicted 

9l)  aU 

by  assuming  -^  = 0. 

The  propellor  seems  to  act  more  as  a "swirl  inducer"  than  a "mixer" 
in  this  near  wake  region.  This  is  supported  by  the  data, in  that  the 
levels  of- the -maximum  temperature  observed  do  not  fall  rapidly  with 
distance  downstream.  If  a severe  mixing  action  were  occuring,  one 
would  expect  a nearly  uniform  temperature  to  appear  at  these  downstream 
locations.  That  this  is  not  the  case  indicates  that  there  is  a more  or 
less  fixed  fluid  element  rotation  occuring.  Future  wake  models  should 
account  for  this  behavior. 

Turb ulence  Obs ervations  in  the  Wake 

Turbulence  measurements  were  performed  in  the  near  wake  at  stations 
Z/D  = 0.33  and  4.0.  At  each  station, axial  and  radial  turbulence  intensi- 
ties, radfaT ’ shear “stre'ss,  and  temperature  fluctuations  were  recorded. 
Figures  No.  (41)  - (47)  are  graphs  of  the  axial  and  radial  turbulence 
intensities  at  these  stations. 

Before  discussing  these  data  a condition  that  existed  which  is 
detrimental  to  the  accurate  measurement  with  a hot-wire  should  be  pointed 
out.  The  injected  air  had  a small  oil  content  which  was  present  because 
of  leaky  compressors.  Examination  of  the  hot-wire  instruments  after  the 
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tests  revealed  that  small  amounts  of  oil  were  deposited  on  the  wires. 

It  is  probable  that  the  electrical  and  thermal  properties  of  the  wires 
were  changed  somewhat  by  this  deposit.  Accordingly,  the  absolute  levels 
of  the  reported  fluctuations  could  be  somewhat  in  error.  Several  measure- 
ments were,  however,  duplicated  at  different  times  and  with  widely  vary- 
ing conditions  such  as  compressor  run-time  (and  thus  oil  vapor  accumula- 
tion) and  the  data  were  seen  to  agree  within  10%.  This  would  indicate 
that  the  error  due  to  the  oil  deposit  is  minimal.  Figures  No.  (40)  - 

(42)  show  that  the  axial  velocity  fluctuations-are  high  (nearly  1 0%.  .oil 

the  free  stream  velocity)  at  the  centernJf  the  wake.  This  is  at  variance 
with  the  data  reported  in  Ref.  No.  (22)  which  show  a minimum  at  this 
point.  The  fluctuations  decrease  rapidly  followed  by  another  peak  which 
occurs  in  the  vicinity  of  the  propel  lor  tips.  This  is  identical  to  the 
trend  of  the  data  in  the  above  reference.  At  Z/D  = 4.0,  Figs.  No. 

(44)  - (46)  show  generally  the  same  distribution  although  the  maximum 
intensities  in  the  wake  center  have  fallen  to  about  half  their  former 
values.  The  axial  turbulence  intensity  in  the  area  of  the  propellor 
tips  for  the  vertical  traverse  of  Fig.  No.  (44)  has  not  decreased  nearly 
as  much.  The  radial  turbulence  intensity  along  all  three  traverses  at 
Z/D  = 0.33  has  the  same  distribution,  with  the  maximums  at  the  center  of 
wake  and  at  the  propellor  tips.  However,  at  Z/D  = 4.0  the  radial  velocity 
fluctuation  has  fallen  to  1 - 2%  of  the  free  stream  velocity  near  the 
tips.  The  fluctuations  along  the  center  of  the  wake  remain  relatively 
high  at  3 - 5%. 

Figures  No.  (43)  and  (47)  are  graphs  of  the  radial  shear  stress 
at  stations  Z/D  = 0.33  and  4.0  respectively.  During  this  point  of  the 
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experimental  program  there  were  equipment  problems  with  the  correlator. 
Consequently,  there  were  several  points  at  which  the  shear  stress  could 
not  be  determined.  The  data  that  we  were  able  to  obtain  indicate 
that  the  highest  values  of  the  shear  stress  occur  in  the  vicinity  of 
high  velocity  gradients,  which  is  to  be  expected.  The  maximum  values  of 
the  shear  stress  do  not  decrease  greatly  from  Z/D  = 0.33  to  4.0. 

Figure  No.  (48)  is  a graph  of  the  temperature  fluctuation  at  Z/D 
= 0.33.  These  results  show  that  the  fluctuations  are  uniformly  very  small 
(<  . 5 % of  Tj  along  any  traverse  until  the  edge  of  the  wake  is  approached. 
The  temperature  fluctuations  in  this  region  are  a maximum  of  approximately 
10.5%  along  the  horizontal  and  e = 45°  traverses  and  a comparatively 
large  14%  along  the  vertical  axis.  Comparison  with  the  data  of  Fig. 

No.  (49)  shows  the  same  behavior.  Maximums  of  V T7^” /T  are  reduced  to 

oo 

less  than  10%  along  the  horizontal  and  e = 45°  traverses  and  to  about 

13%  along  the  vertical  traverse.  Again  the  levels  along  the  vertical 

* 

axis  in  the  vicinity  of  the  propellor  tips  are  larger . The  results 
in  Figs.  No.  (35)  and  (39)  show  that  the  mean  temperatures  are  higher 
in  the  vicinity  of  the  center  of  the  wake  and  much  smaller  at  the  edge. 

The  turbulence  results  indicate  that  the  maximum  temperature  fluctuations 
are  occuring  in  a region  of  relatively  low  mean  temperature  gradients. 

This  is  an  important  result  when  it  is  considered  in  the  light  of  cur- 
rent techniques  used  in  the  theoretical  study  of  compressible 
turbulent  flowfields.  During  the  course  of  such  a study,  at  some 
point  it  becomes  necessary  to  depart  from  the  exactness  of  the  Navier- 
Stokes  equations  and  model  various  terms  containing  the  fluctuating 
components  of  density,  velocity  and  temperature.  An  example  of  such  a 
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term,  for  example,  is  the  turbulent  heat  flux 


qi  s CpP  U!T* 


where  the  subscript  i = 1,  2,  3 denotes  the  coordinate  direction.  The 
usual  practice  is  to  model  the  U|T‘  term  as 

3T 1 
aX ' 


U T 1 

l 


From  this  model,  one  would  expect  a larger  \ T'2  / To  in  the  interior 
3T 

of  the  wake  where  — is  larger  and  comparatively  smaller  fluctuations 
jy 

in, the  outboard  region.  The  experimental  results  of  this  study  do  not 
confirm  such  a model.  There  is,  however,  some  indication  that  has 

dependence  on  -|jp  . The  data  of  Z/D  = 0.33  and  Z/D  = 4.0  ( Figs.  No. 
(48)  and  (49))  show  larger  fluctuations  in  the  tip  area  along  they 
axis  than  in  the  tip  areas  along  the  0 = 45°  traverse  and  X axis.  The 
gradients  in  the  mean  temperature  are  larger  in  the  tip  region  along  the 
Y axis . 


The  questioned?  why  the  temperature  fluctuations  are  concentrated 
only  in  the  tip  regions  is  a difficult  one  to  answer.  On  the  surface, 
it  is  tempting  to  answer  that  the  fluctuations  are  directly  associated 
with  the  vortices  which  are  shed  at  the  tips.  Indeed,  a common  assump- 
tion of  blade  element  theory  to  calculate  propel  lor  performance  is  to 
assume  that  the  circulation  is  constant  along  the  blade  and  consequently 
vortices  are  shed  only  at  the  root  and  the  tips.  Close  examination  of 
Figs.  No.  (48)  and  (49)  show  that  the  temperature  fluctuations  do  tend 
to  rise  at  the  wake  center  in  the  vicinity  of  the  roots  of  the  propellor 
blades.  However,  the  strength  of  circulation  does  vary  along  areal  propellor 
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blade.  In  fact,  it  can  be  proved  that  the  condition  of  constant 
circulation  along  the  blade  length  is  physically  impossible  (e.g.,  see 
Ref.  No.  (24)).  This  means  that  vortices  are  necessarily  shed 
along  the  entire  span  so  the  high  temperature  fluctuations  at  the  tips 
cannot  be  explained  solely  by  vortex  shedding.  The  tip  vortices  will, 
however,  be  of  greater  strength  and  correspondingly  higher  turbulence 
levels  should  be  expected. 

Another  possibility  is  the  effect  of  compressibility.  The  tip 
speeds  associated  with  this  experiment  are  approximately  400  ft/sec 
(M  = 0.35).  Although  this  is  below  M = 0.5,  the  value  of  Mach  number 
above  which  compressibility  effects  are  normally  considered  to  be  import- 
ant, it  is  sufficiently  high  that  some  measurable  compressibility 
effects  should  be  expected.  There  is  the  additional  possibility  that 
the  higher  temperature  fluctuations  at  the  tips  are  the  result  of  eddies 
being  convected  along  the  span  of  the  blade  and  finally  being  shed  at 
the  tips.  This  would  happen  if  the  flow  is  separating  as  it  moves  over 
the  airfoil  surface.  If  this  were  the  case,  a corresponding  loss  of 
thrust  should  be  apparent.  There  was,  however,  no  indication  of  stalled 
propellor  blades  in  any  of  the  flowfield  data. 


V.  CONCLUSIONS  AND  RECOMMENDATIONS 


An  experimental  study  of  the  flow  field  produced  by  a self-propel- 
led, slender  body  of  revolution  in  a temperature  stratified  airstream 
has  been  conducted.  The  experiments  concentrated  mainly  on  the  mean 
and  turbulent  flow  in  the  near-wake  region  although  mean  boundary  layer 
measurements  were  taken  in  order  to  establish  the  initial  conditions 
immediately  upstream  of  the  propel  lor. 

The  mean  flow  results  have  shown  that  the  temperature  stratifi- 
cations, as  produced  in  this  experiment,  has  no  large,  noticeable  effect 
on  the  development  of  the  mean  flow  in  the  near  wake.  The  temperature 
results  show  that  the  fluid  undergoes  a simple  rotation  in  this  near- 
wake reaion  and  that  there  is  no  great  increase  in  mixing  rate  produced 
by  the  propel  lor. 

The  turbulence  measurements  strongly  reinforce  recent  findings 
that  the  highest  levels  of  axial  turbulence  intensity  are  to  be  found 
in  the  immediate  vicinity  of  the  propellor  tips.  Additionally,  the 
current  results  indicate  that  higher  levels  of  turbulence  than  previ- 
ously found  are  present  in  the  wake  center.  These  centerline  fluctu- 
ations appear  to  decay  at  a faster  rate  than  do  the  turbulence 
intensities  at  the  outboard  edge.  They  are,  possibly,  only  a phenomenon 
peculiar  to  the  near-body  region.  The  radial  turbulence  intensities 
are  also  unusually  high  at  the  wake  center.  In  fact,  in  this  case,  the 
magnitudes  along  the  centerline  are  higher  than  those  near  the  propellor 
tips  at  the  most  distant  downstream  station  studied.  The  turbulent 
fluctuations  in  temperature  quite  emphatically  point  out  the  unsteadi- 
ness of  the  flow  in  the  outer  regions  of  the  wake.  For  the  conditions 
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of  this  experiment,  this  was  a vicinity  of  only  small  to  moderate 
gradients  in  mean  temperature.  Temperature  fluctuations  in  the  wake 
in  areas  other  than  the  vicinity  of  the  propellor  tips  are  quite  small. 

This  experimental  program  has  pointed  out  the  relative  ease  with 
which  measurements  of  mean  temperature  and  temperature  fluctuations  can 
be  made.  It  proved  to  be  simpler  to  use  the  straight  hot-wire  in  the 
low  overheat  mode  to  measure  the  temperature  fluctuations  than  to  gather 
the  velocity  fluctuation  data  with  the  cross-wire.  Additionally,  the 
measurement  of  mean  temperature  data  with  a thermocouple  is  very  rapid 
and  uncomplicated.  Accordingly,  it  is  suggested  that  future  experimental 
programs  not  hesitate  to  gather  these  data  in  any  situation  where  temper- 
ature fluctuations  may  be  an  element  of  interest. 

The  results  suggest  some  interesting  points  for  further  study.  The 
influence  of  different  initial  mean  temperature  distributions  should  be 
determined.  The  temperature  fluctuation  distribution  in  the  vicinity  of 
the  propellor  tips  is  quite  pronounced  and  it  would  be  interesting  to 
see  the  effect  of  different  initial  temperature  gradients  and  their 
location  relative  to  the  propellor  geometry.  Second,  the  turbulent 
thermal  boundary  layer  should  be  considered  since  the  incoming  turbu- 
lence levels  should  have  some  effect  on  the  wake  development. 
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Fig.  48  Temperature  Fluctuation  Distribution  at  Z/D  = 0.33 


3.0 


ibution  at  Z/D  = 


APPENDIX  A 
TABULATED  DATA 


3<? 


83 


TABLE  A-l 


Mean  Velocity  and  Temperature  Profiles  on  Body  at  Z'/D  = 2.67 
Injection  Temperature  = 212°F  - 217°F 

Free  Stream  Properties:  T^  72°F,  P^  27.91  in.  Hg,  u^  69.33  ft/sec. 


V ( in. ) 

u/u 

oo 

T/T  (T  = °F) 

oo 

0.0050 

0.510 

1 .458 

0.0194 

0.543 

1 .458 

0.0340 

0.752 

1 .458 

0.0450 

0.763 

1 .458 

0.0586 

0.848 

1 .458 

0.0722 

0.955 

1 .458 

0.0858 

1 .000 

1 .458 

0.0994 

1 .010 

1 .458 

0.1260 

1 .022 

1 .458 

0.1540 

1 .023 

1 .444 

0.1820 

1 .025 

1 .444 

0.2360 

1 .026 

1 .431 

0.2960 

1 .026 

1 .417 

0.3560 

1 .025 

1 .375 

0.4160 

1 .023 

1 .347 

0.4760 

1 .021 

1 .319 

0.5660 

1 .019 

1 .278 

0.6560 

1 .015 

1 .222 

0.7460 

1 .014 

1 .20? 

0.8360 

1 .012 

1 .181 

0.9560 

1 .009 

1 .125 

1 .0760 

1 .006 

1 .097 

1 .1960 

1 .005 

1 .069 

2.0800 

1 .000 

1 .000 
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TABLE  A-2 


Mean  Velocity  and  Temperature  Profiles  on  Body  at  Z ‘ /D  = 4.67 
Injection  Temperature  = 212°F  - 217°F 

Free  Stream  Properties:  T = 69°F,  Pc  =27.81,  u = 69.26  ft/sec 

oo  ^co  oo 


V (In.) 

u/u 

oo 

T/T  (T  = °F ) 

oo 

0.0050 

0.476 

1 .435 

0.0194 

0.500 

1 .435 

0.0340 

0.707 

1 .435 

0.0450 

0.698 

1 .435 

0.0586 

0.764 

1 .435 

0.0722 

0.790 

1 .435 

0.0858 

0.808 

1 .435 

0.1260 

0.864 

1 .435 

0.1540 

0.876 

1 .435 

0.1820 

0.922 

1 .435 

0.2360 

0.965 

1 .420 

0.2960 

1 .007 

1 .406 

0.3560 

1 .025 

1 .391 

0.4160 

1 .028 

1 .377 

0.4760 

1 .026 

1 .348 

0.5660 

1 .024 

1 .319 

0.6560 

1 .022 

1 .290 

0.7460 

1 .019 

1 .246 

0.8360 

1 .017 

1 .217 

0.9560 

1 .015 

1 .174 

1 .0760 

1 .012 

1 .130 

1 .1960 

1 .009 

1 .087 

2.1100 

1 .000 

1 .000 
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TABLE  A-3 


Mean  Velocity  and  Temperature  Profiles  on  Body  at  Z'/D  = 6.67 

Injection  Temperature  « 212°F  - 217°F 

Free  Stream  Properties:  T = 69°F,  Pc  = 27.91,  u = 69.14 

OO  jC®  oo 


Y'(in.) 

u/u 

oo 

T/T  (T  = 

oo 

0.0050 

0.440 

1 .406 

0.0194 

0.570 

1 .406 

0.0340 

0.638 

1 .406 

0.0450 

0.650 

1 .406 

o.„::: 

0.699 

1 .406 

0.0722 

0.704 

1 .406 

0.0858 

0.740 

1 .406 

0.1260 

0.777 

1 .406 

0.1540 

0.801 

1 .406 

0.1820 

0.838 

1 .391 

0.2360 

0.873 

1 .391 

0.2960 

0.927 

1 .391 

0.3560 

1 .954 

1 .391 

0.4160 

0.979 

1 .377 

0.4760 

1 .999 

1 .362 

0.5660 

1 .009 

1 .333 

0.6560 

1 .012 

1 .290 

0.7460 

1 .013 

1 .275 

0.8360 

1 .014 

1 .246 

0.9560 

1 .013 

1 .203 

1 .0760 

1 .010 

1 .159 

1 .1960 

1 .008 

1 .130 

2.0900 

1 .000 

1 .000 
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TABLE  A-4 


Mean  Velocity  and  Temperature  Profiles  on  Body  at  Z'/D  = 8.67 
Injection  Temperature  = 212°F  - 217°F 

Free  Stream  Properties:  T = 68°F,  Pc  = 27.90,  u = 69.08  ft/sec 

CO  J QO  ro 


V'( in.) 

u/u 

oo 

T/T  (T  = 

oo 

0.0050 

0.400 

1 .397 

0.0194 

0.450 

1 .397 

0.0340 

0.545 

1 .397 

0.0450 

0.620 

1 .397 

0.0586 

0.680 

1 .397 

0.0722 

0.700 

1 .397 

0.0858 

0.720 

1 .397 

0.1260 

1 .751 

1 .397 

0.1540 

1 .752 

1 .382 

0.1820 

0.795 

1 .382 

0.2360 

0.826 

1 .368 

0.2960 

0.871 

1 .368 

0.3560 

1 .905 

1 .368 

0.4160 

1 .926 

1 .353 

0.4760 

0.952 

1 .338 

0.5660 

0.976 

1 .294 

0.6560 

0.990 

1 .265 

0.7460 

0.996 

1 .250 

0.8360 

1 .999 

1 .221 

0.9560 

1 .004 

1 .176 

1 .0760 

1 .001 

1 .147 

1 .1960 

1 .000 

1 .118 

2.3600 

1 .000 

1 .000 
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TABLE  A-5 


Mean  Velocity  and  Temperature  Profiles  on  Body  at  Z'/D  = 10.0 
Injection  Temperature  = 212°F  - 217°F 

Free  Stream  Properties:  T = 68°F,  Pc  = 27.90,  u = 69.08  ft/sec 

co  joo  oo 


Y * (in.) 

u/u 

oo 

— 1 

8 

— i 

0.0050 

0.380 

1 .368 

0.0194 

1 .450 

1 .368 

0.0340 

0.500 

1 .368 

0.0450 

0.540 

1 .368 

0.0586 

0.617 

1 .368 

0.0722 

0.630 

1 .368 

0.0858 

0.647 

1 .368 

0.1260 

1 .678 

1 .353 

0.1540 

1 .740 

1 .353 

0.1820 

1 .770 

1 .353 

0.2360 

0.816 

1 .353 

0.2960 

0.827 

1 .353 

0.3560 

1 .846 

1 .353 

0.4160 

0.873 

1 .338 

0.4760 

1 .906 

1 .324 

0.5660 

0.940 

1 .309 

0.6560 

0.970 

1 .294 

0.7460 

1 .997 

1 .265 

0.8360 

1 .012 

1 .235 

0.9560 

1 .013 

1 .206 

1 .0760 

1 .010 

1 .162 

1 .1960 

1 .009 

1 .147 

2.4000 

1 .000 

1 .000 
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TABLE  A- 1 3 


Temperature  Fluctuations  in  the  Wake  at  Z/D  =0.33 

Injection  Temperature  = 172°F  - 178°F,  Q = 1 .0  in.  H„0,  T = 41°F 

L 00 


X(  in . ) 


0.00 
0.00 
0.00 
0.00 
0.00 
0.00 
0.00 
0.00 
0.00 
0.00 
0.00 
0.00 
0.00 
0.00 
0.00 
0.00 
0.00 
0.00 
0.00 
0.00 
0.00 
0.00 
0.00 
0.00 
-1.50 
-1.00 
-0.85 
-0.50 
-0.25 
0.00 
0.25 
0.50 
0.75 

1 .00 

1.50 

2.00 

2.50 
2.75 

3.00 
3.25 


Y(in.) 

PSoo(in.) 

V T‘2 

T'oo 

-1.50 

27.57 

0.038 

-1 .00 

27.57 

0.031 

-0.75 

27.57 

0.022 

-0.50 

27.57 

0.011 

-0.25 

27.57 

0.014 

0.00 

27.57 

0.017 

0.25 

27.58 

0.026 

0.50 

27.58 

0.032 

0.75 

27.58 

0.024 

1 .00 

27.58 

0.030 

1 .50 

27.58 

0.047 

2.00 

27.58 

0.081 

2.25 

27.53 

0.110 

2.50 

27.59 

0.122 

2.75 

27.59 

0.130 

3.00 

27.59 

0.134 

3.25 

27.59 

0.138 

3.50 

27.59 

0.128 

3.75 

27.59 

0.116 

4.00 

27.59 

0.083 

4.50 

27.58 

0.030 

5.00 

27.58 

0.004 

5.50 

27.58 

0.001 

6.00 

27.58 

0.001 

0.00 

27.57 

0.032 

0.00 

27.57 

0.021 

0.00 

27.57 

0.019 

0.00 

27.57 

0.025 

0.00 

27.57 

0.025 

0.00 

27.57 

0.018 

0.00 

27.57 

0.020 

0.00 

27.57 

0.018 

0.00 

27.57 

0.020 

0.00 

27.57 

0.026 

0.00 

27.57 

0.043 

0.00 

27.57 

0.076 

0.00 

27.57 

0.091 

0.00 

27.57 

0.100 

0.00 

27.57 

0.103 

0.00 

27.58 

0.096 
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TABLE  A-13  (Continued) 


X (in. ) 

Y(  in. ) 

PSoo(in.) 

vv7 

T 

oo 

3.50 

0.00 

27.58 

0.089 

3.75 

0.00 

27.58 

0.085 

4.00 

0.00 

27.58 

0.066 

4.25 

0.00 

27.58 

0.057 

4.50 

0.00 

27.58 

0.048 

5.00 

0.00 

27.58 

0.053 

5.50 

0.00 

27.58 

0.044 

-1.50 

-1.50 

27.58 

0.053 

-1.00 

-1.00 

27.58 

0.044 

-0.50 

-0.50 

27.58 

0.031 

0.00 

0.00 

27.58 

0.021 

0.50 

0.50 

27.58 

0.031 

1 .00 

1 .00 

27.58 

0.061 

1.50 

1.50 

27.58 

0.116 

1 .75 

1 .75 

27.58 

0.124 

2.00 

2.00 

27.58 

0.109 

2.25 

2.25 

27.58 

0.095 

2.50 

2.50 

27.58 

0.071 

3.00 

3.00 

27.58 

0.013 

3.50 

3.50 

27.58 

0.001 
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TABLE  A- 1 4 


Temperature  Fluctuations  in  the  Wake  at  Z/D  = 4.0 

Injection  Temperature  = 169°F  - 174°F,  Q = 1 .0  in.  Ho0,  T = 41°F 

L 00 


X( in. ) 

Y(in.) 

PsJin.) 

T 

00 

0.00 

-1.50 

27.54 

0.024 

0.00 

-1.00 

27.54 

0.022 

0.00 

-0.50 

27.55 

0.019 

0.00 

-0.25 

27.54 

0.023 

0.00 

0.00 

27.54 

0.028 

0.00 

0.25 

27.55 

0.032 

0.00 

0.50 

27.55 

0.031 

0.00 

0.75 

27.55 

0.028 

0.00 

1 .00 

27.55 

0.027 

0.00 

1.25 

27.55 

0.028 

0.00 

1 .50 

27.55 

0.029 

0.00 

1 .75 

27.54 

0.034 

0.00 

2.00 

27.54 

0.041 

0.00 

2.25 

27.55 

0.065 

0.00 

2.50 

27.56 

0.086 

0.00 

2.75 

27.55 

0.101 

0.00 

3.00 

27.55 

0.123 

0.00 

3.50 

27.56 

0.136 

0.00 

3.75 

27.55 

0.121 

0.00 

4.00 

27.55 

0.108 

0.00 

4.50 

27.55 

0.059 

0.00 

4.75 

27.56 

0.024 

0.00 

5.00 

27.56 

0.002 

0.00 

6.00 

27.56 

0.001 

-1.50 

0.00 

27.57 

0.050 

-1.00 

0.00 

27.57 

0.044 

-0.50 

0.00 

27.57 

0.029 

-0.25 

0.00 

27.57 

0.025 

0.00 

0.00 

27.57 

0.025 

0.25 

0.00 

27.57 

0.025 

0.50 

0.00 

27.57 

0.026 

0.75 

0.00 

27.57 

0.032 

1.25 

0.00 

27.57 

0.040 

1.50 

0.00 

27.57 

0.027 

1.75 

0.00 

27.57 

0.055 

2.00 

0.00 

27.57 

0.064 

2.50 

0.00 

27.57 

0.079 

3.00 

0.00 

27.58 

0.086 

3.25 

0.00 

27.58 

0.088 

3.50 

0.00 

27.58 

0.072 

L_ 
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TABLE  A- 1 4 

(Continued) 

X (in. ) 

Y(in.) 

Pc  (in.) 

b°° 

Vti2 

Too 

4.00 

0.00 

27.58 

0.042 

4.50 

0.00 

27.58 

0.022 

5.00 

0.00 

27.58 

0.028 

-1.50 

-1.50 

27.58 

0.050 

-1.00 

-1.00 

27.58 

0.042 

-0.50 

-0.50 

27.58 

0.020 

0.00 

0.00 

27.58 

0.026 

0.50 

0.50 

27.58 

0.026 

1 .00 

1 .00 

27.58 

0.033 

1.50 

AM]  .50 

27.58 

0.081 

1.75 

1.75 

27.57 

0.094 

2.00 

2.00 

27.57 

0.086 

2.50 

2.50 

27.57 

0.037 

2.75 

2.75 

27.57 

0.022 

3.00 

3.00 

27.57 

0.001 

